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Abstract 

We obtain the Ward identities and the gauge-dependence of Green's functions in non-Abelian gauge 
theories by using only the canonical commutation relations and the equations of motion for the Heisen- 
berg operators. The consideration is applicable to theories both with and without spontaneous sym- 
metry breaking. We present a definition of a generalized statistical average which ensures that the 
Fourier images of temperature Green's functions of the Fermionic fields have only even-valued frequen- 
cies. This makes it possible to set up a procedure of gauge-invariant statistical averaging in terms of 
the Hamiltonian and the field operators. 

1 Introduction 

The study of the effects of spontaneously-broken gauge theories in statistical physics [T] has raised the 
problem of finding a proof of the gauge- invariance of physical results in gauge statistical physics. Formally, 
such a proof can be carried out by analogy with quantum field theory, by using a representation of statistical 
averaging with the help of functional integration, which has been presented in [2]. However, there exist 
some specific calculations of physical effects that have an appearance of being gauge-dependent. This 
circumstance may cast a shadow on the applicability of the functional approach to statistical physics. In 
particular, one may raise the question as to the validity of a non-local change of variables in the functional 
integral for the partition function, which has to be made in the course of the usual proof of gauge- invariance, 
as well as in the process of deriving the Ward identities. 

It appears useful, therefore, to deduce the Ward identities and the gauge-invariance of physical results 
in the framework of the operator formalism of quantum field theory, i.e., by using only the Heisenberg 
equations of motion and canonical commutation relations. This has been done in the present article for 
both field theory (Sections 3, 4) and statistical physics (Section 5). 

It is essential that the construction of a theory requires to sum up a gauge-invariant Lagrangian, 
Lq, not only with a gauge-fixing term, but also with an additional Lagrangian, Lc, responsible for an 
interaction of fictitious particles with the gauge field. For those gauges that are usually applied in quantum 
electrodynamics, the fictitious particle is free, so that the theory can be set up without the term Lc- 
In the case of non-Abelian gauge theories, it is well-known that there arise some additional diagrams 
(with respect to the Feynman diagrams) that effectively describe an interaction of the gauge field with 
the fictitious particles. We choose the Lagrangian Lc in such a way that, on the one hand, it implies the 
necessary additional diagrams, and, on the other hand, it admits a canonical quantization. 

The deduction of the Ward identities and the proof of gauge-invariance that are based on the equations 
of motion for the Heisenberg fields can also be useful for other purposes, such as a description of gauge 
theories in the framework of Zimmermann's normal product ,4 . 

The fictitious particles described by the Lagrangian Lq are scalars, but, at the same time, they are 
Fermions. Thus, the usual definition of statistical average leads to such a Fourier-image of the temperature 
Green function of these particles that contains only odd- valued frequencies. However, as will be shown 
in Section 5, the gauge-invariance of the partition function demands that the Green function of fictitious 
particles should contain only even- valued frequencies (whereas the Green function of Bose particles should 
contain odd- valued frequencies). This allows one to pose statistics in terms of an operator formalism by 
using the Hamiltonian and the field operators. Physical quantities, and, in particular, the partition function. 
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prove to be gauge- invariant. In the case of gauges that ehminate the fictitious particles, the generahzed 
definition of the partition function becomes identical with the usual definition. 

The present consideration has been made for a sufficiently large class of gauge conditions and is appli- 
cable to theories both with and without spontaneous symmetry breaking. 



2 General formulas 

In this section, we present some general formulas that we use in the following sections. 

We examine a Lagrangian of a renormalizable theory of a gauge-invariant interaction of spinless and 
Fermionic fields with a gauge field of a general kind: 

L = Lo + ^t''aabt' + Lc, (2.1) 

+ ^7''0^m + .9^°^m)'^-^('/^'^) ' (2-2) 

where (pt is a multiplet of Hermitian spinless fields; Vi is a multiplet of Fermionic fields; V {tp, ip) is an 
arbitrary Lagrangian of interaction of the fields ip, f that obeys the requirement of gauge-invariance; GJ^^j^ 
is the strength tensor of a Yang-Mills field A^: 

= d^At ~ a.^;: + gr'^AlAl , (2.3) 

while /"^"^ are structure constants of an invariance group G; the matrices F" and r" are the generators of 
the transformation group for the fields cpi and ipi, with the commutation relations 

[F^F**] = (2.4) 

and obey similar relations. The matrices F° are antisymmetric and have purely imaginary values. 

The last two terms in (|2.ip describe subsidiary conditions that are necessary for the possibility of a 
canonical description of the theory (as well as for the existence of a perturbation theory). Let us choose 
the functions t"" in the form 

r = ^^"d^A'^, + , (2.5) 

where >f^'' — x'^^ is a matrix that has numerical elements; the matrices will be described below. 

The Lagrangian Lc describes an interaction of fictitious particles (Fermionic scalars), C°, C"*"", with 
the fields (pi and implies the well-known additional diagrams in quantum theory [5^, namely, 

Lc = -d^C+'^x^'^Vfc'' + igC+''x^T'l^(pjC' , (2.6) 
^ d^S'''' + gr^'^Al . 

The matrix xf is chosen as follows: let the vacuum mean value of the field ipi be non- vanishing, 

(0|^,|0)=e.- (2.7) 
The set of generators Tf^ can always be spht in two groups, F° = (F™, F'), such that 

F'^==0, F"C = iC7^0, (2.8) 
where the vectors ^"^ are orthogonal. Then xf must obey the conditions 

>cU^ = 0■ (2.9) 

Besides, either (for a fixed a) 

xf=0, (2.10) 

or 

xtCr^O, (2.11) 

with a certain m. In other respects, the matrix xf is arbitrary. We do not discuss any restrictions that may 
be imposed on xf by conditions (|2.9p - (|2.1ip . Notice that conditions (|2.9p and (|2.1ip are necessary only 
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for the possibility of passing to unitary gauges of the kind x°^Lpi — 0. By themselves, they are unnecessary 
to provide the possibility of a canonical description. 

In principle, one can choose a subsidiary condition of a yet more general kind, being compatible with 
renormalizability: 

t ^ dA + dip + A + + + Aip + ip^ . 

This, however, will only imply some obvious modifications of the reasonings to be presented below. 

Since we examine a theory which admits a possibility of spontaneous symmetry breaking, let us now 
introduce some spinless fields, whose vacuum mean value is zero: 

P^ = U + <^^. (0|a,|0> =0. (2.12) 

The Lagrangian Lq is invariant with respect to gauge transformations whose infinitesimal form is given by 

A1 ^ A'^ = Al + VfA" , V;:" = 8,5'^' + gr'A^^ , (2.13) 
^ V' = V' + igA^T^^P , tp' ^tp + igA^Vip , (2.14) 

where A° are infinitesimal parameters of the gauge transformations; they depend on the coordinates. For 
the field a, transformation (|2.14p reads as follows: 

(T ^ o-' = cr + i^A^r'^^ + igA'^rV , (2.15) 

i.e., a transforms in a non- homogenous manner. The invariance of Lq with respect to (|2.13p - (l2.15|) implies 
the identities 

_ v""!^ + Z9^r>, - zg^r^V, - ^9i^.r^/-^ ^ . (2.16) 

Let us now proceed to the canonical quantization of the theory. The momenta of the fields AJ^, (p and 
■0 are defined as usual: 

TT^ : n^^o = aabx"°t'' , n"-'' = -0"^°'' + Uab^^^t^ , (2.17) 
: = a, ~ igVlf^A^pk , (2.18) 
: n, = (2.19) 

Using these relations, one can obtain the derivatives of the fields with respect to time: 

1 



t = — rjra TT ' ,a Offoc = 0^ , (2.20) 



A>0 
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Ok 

i'^^fe ^ —n'''° - tt"'*^ + V''''^'^A''^° , (2.21) 

^0 = ^ - ^""'^'^M - + ^o.tt''-^-- (2.22) 

OA; 

a = Tl, + igri^Alp,. (2.23) 

The canonical variables obey the usual equal-time commutation relations 

[Al, ^"■'^j - iS'^H^, , [a,,U,] = i5,, , [0„ Hj] = i5,j . (2.24) 

To find the canonical momenta conjugate to the fictitious fields, as well as the corresponding anticom- 
mutation relations, one can use Schwinger's action principle [B]; see Appendix A, where it has also been 
demonstrated that these expressions are, in fact, the only possible ones. As a result, we find 

Tl^'c ^ -><°'"a^C+" , = x^f^'^fC^ , (2.25) 

= ^ (n^ + >r°'^9feC+") , 

C"* ^ _L (n;^.+ - gx^^p^'^AlC'^ - x°^\/fC'') , (2.26) 

{c^n^} = {c+",n^+} ^ifj"''. (2.27) 
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The other anticommutators are equal to zero. The corresponding Hamihonian reads 

H = n'^'^'A^^ + n.a, + + n^c" + n;5.+c'+" -l = 



2(^00)^ V / 2 

1 „„ 1 



Pi" 



+ (^^^^■'' - ^^^J^) Vf - i5C+"<r,^(^jC'' . (2.28) 

One can easily see that the canonical equations that follow from Hamiltonian (|2.28p are identical with the 
Lagrangian equations. For the fictitious particles, this has been verified in Appendix A. 

We need an expression for tt"'" in terms of the canonical coordinates and momenta. It can be found 
with the help of the following relation, that holds true for any operator Q: 

Q=^i[H,Q] . (2.29) 

We have 

•a,0 _ ^ Oipf afl ^ Oiyjab b,0 , -ryab b,k _ 

- igU.Tt^ip, + #70rV + gH'cf'^'c' . (2.30) 

From (|2.16p , which holds identically, it follows that on the equations of motion relation (|2.16p is valid also 
for the quantity L — Lq, that is, 

= -gnt'^d^ (9^C+V'""^C''^) " g"^ Af x^"" d^C+'^ J'^''^ + 

+ g^C+''x\T%T';^^kC'' , Af - r^A^ , (2.31) 
where T"-^ stands for the operator 

T'^fc = a^^^A-^vf + igx^T^Vj ■ (2.32) 
The equations of motion for the fictitious fields read as follows: 

rpabfjb ^ (j+bj,ba _ (-2 

3 Ward identities 

In this section, we deduce the Ward identities by using only the equations of motion and canonical com- 
mutation relations. 

Consider the vacuum mean value 

Zc^c^. = (0 \TC'' (y) {x) exp {iQ)\ O) , (3.1) 

where 

g = J duQ{u)^ j du [J^^'^Al + J.ipi +r]ij+ Iprj) . (3.2) 
An arbitrary operator P (z) obeys the relation 

9^'^ {P (2)> ^(P{z))+i j du5 (20 - u^) {[P {z) , Q (u)]) + 
+ (0 \T {exp {iQ) 5 {zo - yo) [P (z) , C (y)] C+'' {x) + 

+ {y) 6 {zo - xo) [P {z) , C+« {x)] } | O) , (3.3) 
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where the following notation has been used: 

{P [z)) ^ {0\TP [z] [y) C+'^ (x) exp {^Q) |0) . (3.4) 
Because of the fact that all the fields (anti) commute at equal times, we have 

t^' (z) Zcc+ ^ (tUz)) . (3.5) 
Expressions of the form P (A, ip, -0, ip) Zq(j+ imply that the function P is subject to the replacement 



. 1 ^ 
A -— , etc. 

I oJ 



In particular, 



tf (z) Zcc^ = x^'-S^^) {Af (z)) + ^{ (z)) . 
Further, relations ((2l7|) - ((2:271) . Ip^O)) . (jSlS)) lead to 

(z) Zcc+ = {i^ (z)) + 4o (^) ^cc^ ' 

f (z) Zcc+ = (t^ iz)) + ia^^j^^o (z) Zee. - -^o.f (Q^ (z)) 
4o ^^^^^.Z'-^o (z) + 4o >^°'V^ J=-° (z) + VgV-° (z) 



(3.6) 
(3.7) 

(3.8) 



MO' 



5iz-y)-^af'f''''Zc.'c 



^00" •' -^C'''C+'- ■ 

In p.9p . the expression Qji denotes an infinitesimal variation of the term with the sources, 

Qr (z) = -Vl'r-^ (z) + igJ, (z) r^^.^, (z) + igrJiz) r'^p (z) - ig^ (z) t\ (z) . 

Relations p.Sp and p.9p allow one to obtain the following: 

(z) ay,T^^Zcc+ = (i^ (^) a/cT^") - 
-(Q$,^(z))+^<5(z-y)/'"''^'Zc.,c:+- 



(3.9) 



(3.10) 



(3.11) 



We now use relation (|2.31l) in the first summand of the r.h.s. of (13. lip , which is a valid operation, since all 
the derivatives are already under the symbol of T-product; we assume b = d, y = z and take an integral 
over y, as well as a sum over d, namely. 



I dytf{y)af,T^''Zcc+=- J dy {Q%Xy)) + 



dy 



-g^i'^'d^ {0.0+" (y) f'^Cf (y)) - g^Af (y) ^.^^•'d.C+f (y) /^^"C" (y) + 



g^C+\y)^tvi^%^, {y)Ci (y) 



(3.12) 



Using the equations of motion for C and C"*" (|2.33p . as well as the anticommutativity of the operators C, 
one can prove (see Appendix B) that the second summand in the r.h.s. of (j3.12p is equal to zero. Then 
(|3.12p takes the form 

(3.13) 



j dyt' (y) abcT-^Zcc^ = - (^J dyQi (y) 



In view of expression (|2.32p , one can easily see that the derivatives in all the terms T'^'^ commute with the 
symbol of T-product, with the exception of the term do>c°'^\7^. By virtue of (|2.33p . (|2.25p . (|2.27p . we find 



T^'^Zcc+ = iS^^6 {y-z)Z., 
where Z stands for the vacuum mean value 

Z= (0|Texp(iQ)|0) . 



(3.14) 



(3.15) 
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As a result, the Ward identity for the function Z takes the form 



o-aht (x) = iiO 



T dyQ'kiy)C'>iy)C^^x) 



In order to present p.l6p in the usual form, we utilize a relation that follows from p.l4p . 

P {A, if, ^) Zcc+ = iD'''' {y, x) P (A, ^,^)Z, T^'d"'' = 5^ , 



(3.16) 



(3.17) 



for an arbitrary function P. This relation allows one to present (|3.16p in the usual form which can be found 
in the literature: 

aaht''{=^)+ f dyQ''Riy)D''^iy,x) Z = 0. (3.18) 



For the first time, an identity of the form p.lSp has been obtained by Fradkin [7] for an Abelian theory, 
as well as by Slavnov [8J and Taylor 9 for a non- Abelian gauge theory. 

Identity (j3.16p has a simple meaning. The substitution of a new field t (x) into this Green function is 
equivalent to the sum (over the number of fields in the initial Green function) of Green's functions that do 
not contain the field t (x) and are deduced from the initial Green function with the help of an infinitesimal 
gauge transformation ()2.13p ~ ()2.15p of one of the fields with the gauge function A ^ CC+. In case one 
of the fields in the Green function is a gauge-invariant operator, an insertion of the field t (x) into such a 
function yields zero. 



4 Gauge-dependence of Green's functions 

In this section, we find a relation between Green's functions in the gauges Uat, x^'^ , >ff and 

Ciab = aab + 5aab , ^^"^ = x''" + S x''" ,J(f = xf+ (5< . 

As a preliminary step, we make the following remark. In the framework of canonical formalism, the 
generating functionals (|3.ip . (|3.15p . or, equivalently, the Green function, are computed by making use of the 
vertices determined by i?int and also with the help of non-covariant propagators, being the actual vacuum 
mean values of the T-products of free fields. 

In case the Hamiltonian is quadratic in its momenta, one can pass (due to Wick) to an effective dia- 
grammatic technique!! in which the vertices are determined by — iint; the field propagators are determined 
by A~^ (A being a differential operator that enters the free Lagrangian L ^ ^ipAtp), and, besides, there 
may also appear some additional vertices. 

The additional vertices are determined by the matrix present in those summands of the Hamiltonian 
that are quadratic in momenta ^ . In the case under consideration, the effective Lagrangian that determines 
the additional diagrams equals to 

-^(5(0)Splnaa6, (4.1) 

i.e., there are no additional diagrams. Thus, the generating functionals (j3.ip . (|3.15p can be presented in 
the form 

Z = exp Q,5 (0) Sp In Zw , (4.2) 

Zcc+ = exp Q(5 (0) Sp Inaab^ Zcc+w , (4-3) 

where Zw and Zcc+w are now computed with the help of the usual Feynman diagrams. 

It is easy to demonstrate (Appendix C) that the variation of Green's functions with respect to the 
parameters entering the Lagrangian L is given by an insertion of the "field" i J dx SL into Green's functions. 
Due to the fact that the vacuum mean value ^ depends on the gauge parameters, we have 

^Z = i J dx SLZ + i J dx6^i(^^ + Z+^S{Q) a''''5aabZ , (4.4) 



^One can pass to this diagrammatic technique also in the general case. However, in case the Hamiltonian is more than 
quadratic in its momenta, one cannot find the additional diagrams in a manifest form. 
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where 5L denotes a variation of the Lagrangian with respect to those gauge parameters that enter the 
Lagrangian in a manifest way: 



dx 6L = i / dx 



1 



e [ ^a'"'Sacdt'^ + St'' ] aab + C+'^(5T°^C'' 



(4.5) 



The second term in (|4.4p equals to zero, because it represents one of the equations satisfied by Z as a 
function of the sources. This fact foUows from the Euler equations for the fields, as well as from definition 
p.l5p and from the canonical commutation relations. 
Let us now use the Ward identity in the form (|3.16p : 



SZ = - I dx { -a'^'Scbt' (x) + St" (x) 



Tcxp(*Q) / dyQj,{y)C''{y)C+-{x) 



+ ^( 



Tcxp(iQ) / rfx (a;) r"''C'' (x) 



0) + -,5(0) a'^'SaabZ. 



(4.6) 



With allowance for definitions (|2.5p , (|2.32p and (|3.10p , the sum of the second term and the part of the first 
term which contains St in the r.h.s. of (|4.6p yields 



Texp (iQ) I dxdy Q^j^ (x) C° (x) C+'' [y) St" (y) 



(4.7) 



Let us recall that we compute Z with the help of Feynman diagrams, so that the time derivatives of t, 
St, ST and Qn should be regarded as commuting with the symbol of T-product. However, relation (|3.17p 
remains valid. Due to this fact, we finally have 



SZ 



dxdyQ'kix)D'"'ix,y)A (y) Z , 



A" (x) ^ -a^'Habct'' {x) + Sf [x) . 



(4.8) 
(4.9) 



Relation (j4.8p has a simple meaning as well: an infinitesimal change of the gauge parameters in the 
Green function is equivalent to an infinitesimal gauge transformation of each of the fields with the gauge 
parameter A ~ DA. In particular, (non-renormalized) Green's functions of gauge-invariant operators are 
gauge-independent . 

Concluding this section, notice that one could remain in the framework of the Dyson T-product and 
use the variation SZ (Appendix C) in the form 



SZ ■ 



-i ( 



T j dxSHexp {iQ) 



This would only lead to some more tedious calculations, that would naturally leave relation (|4.8p unaltered 



5 Gauge invar iance of partition function 

In this section, we examine the problem of defining the partition function in gauge theories. 
Notice that if one defines the partition function as usual, 

Z = Spe-(^" (5.1) 

[where Sp is defined (in a space with indefinite metric) so that a cyclic permutation of operators under the 
symbol of Sp is admissible], then it is not gauge- invariant, since this case does not comply with the Ward 
identities p.l6p . which provide a basis for the gauge-invariance of physical quantities. 
Indeed, in the simplest case of free electromagnetic field, described by the Lagrangian 

L = -\fI, - d,C+d^C+^ {d.A^f , (5.2) 
identity (|3.16p for the temperature propagator of the field takes the form 

(A^A,) = ~d, {CC+) , (5.3) 
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where we have introduced the notation 



(Q) ^Spe-^'^TrQ (5.4) 

for any operator depending on the "temperature time" r [lOj . The right- and left-hand sides of relation 
(|5.3p can be computed directly, and, in both cases, they are given by the kernel d,^-^. Nevertheless, the 
Fourier-image of the l.h.s is known to contain only even-valued frequencies, whereas the Fourier-image of 
the r.h.s., just as Green's functions of any Fermionic operators, contains only odd-valued frequencies (see 
also Appendix D). Therefore, relation (|5.3p does not take place. 

Since the operator structure and diagrammatic technique for statistical temperature Green's functions 
are known [10] to be completely analogous to the corresponding expressions of field theory, the derivation 
of the Ward identities and gauge properties for temperature functions should be carried out in complete 
analogy with the corresponding calculations of Section 4. 

There is, however, an operation that should be examined in more detail. This is integration by parts, 
which has been used several times in Section 4 (see below). Integration by parts is valid in case the 
corresponding vertex conserves momentum (i.e., the sum of the momenta of all the fields in a given vertex 
equals to zero). Within the temperature techniques, this condition holds true only in case the sum of 
frequencies of all the fields in a vertex is even-valued [10]. As regards the problem we discuss here, the 
"criminal" cases of integration by parts arise as one passes from (|3.12|) to (|3.13p . that is, as one shifts the 
action of a derivative from the field t to the field C in the l.h.s. of (|3.12p . and also as one proves the fact 
that the second term in the r.h.s. of (|3.12p is equal to zero. Given this, one encounters vertices of the kind 

J dydtC, J dyC+CC, j dyC+CA, etc. (5.5) 

One can observe that integration by parts is valid in case the fictitious Fermionic field C (as well as 
the field A ) contains only even-valued frequencies. 

Thus, the partition function defined by formula (j5.4l) is not gauge-invariant. In order to provide gauge- 
invariance, it is also necessary that the Green function of a fictitious Fermionic field should contain, never- 
theless, only even-valued frequencies. Besides, the Wick theorem must also be valid, thus making it possible 
to deduce relation (|4.4p : see also Appendix C. Each of these conditions can be fulfilled. 

Let us recall that the operator 

S = e-'^^+'^'^'^'T, exp (Q) , (5.6) 
dT / d^x ( J'^''' -I- J,ip, + riij + V??) (5.7) 
can be presented in the form 



S = e-^^+'^^'^'T.exp |^g(") - J (^il?t)) • 



(5.8) 



The dependence of the Heisenberg operators in (|5.6p on the temperature parameter r is given by the 
equation 

^A1=[{H-tJi.,N,),Al\ . (5.9) 

For any other operators, it is determined in a similar way and is given by a formal replacement it t. In 
(|5.8p . each of the operators is free, and its dependence on r is determined by the equation 

^A';,= [{Ho-fi.N,),A;] . (5.10) 

We have also introduced some terms with chemical potentials (in case they are necessary), assuming that 
Ni are gauge-invariant. 

Consider some matrix element S (which will be denoted by ( )). It can be presented in the form 

(S) = eXp(-i7int) {S)o\ja^J^^ga^g + a^ga^g+a^o ' (^-^l) 

In (|5.1ip . it is implied that the operator acting on {S)q is subject to the replacement 

A^-^,etc, (5.12) 
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and the function (5)q is defined as follows: 

(5>o = (e-'^^-'+'^^'^-r.exp (q(") + q'"^)) , (5.13) 

J^^ dT J (fx ill'K''^^' + J.n, + e+^-c + c+^-e'' + e'^°'n}.+ + n^^) . (5.14) 

All the operators in (|5.f 3p are free. The matrix elements (|5.1ip . (|5.13p obey the Wick theorem provided 
that the free Green functions determined by (|5.13p [i.e., the coefficients of the series expansion of (|5.13p in 
the powers of sources] should decompose into a product of two-point Green's functions. In other words, 
expression (|5.13p as a function of sources should be given by an exponential of a quadratic form with respect 
to the sources. The following matrix element possesses the required property: 

(^)o = En^r (n|e-'^^o+/5--^-T.exp (q(°) +Q'"')|n) , (5.15) 

n i 

where 

|n) = |ni) ® Inz) (g) . . . . (5.16) 

\ni) are n-particle states, normalized by ±1, of the particles of i-ih type, being eigenstates of the free 
Hamiltonian Hq and of the Hamiltonian Hq — fiiNi; the numbers Xi depend on the type of a particle (in 
principle, they may also depend on the momentum of a particle). In Appendix D, it is shown that (|5.15p 
equals to 



(S), = Z^;-^ exp [^-J^D^pJ^j exp (6 J.) , (5.17) 
where 

4'-' = (^)olj==o ' (5-18) 

while J" denotes the set of all the sources: 

J" = {j; , r-'^ ,J,,J,,rj,rjJ\d-J^\ 9+^} , (5.19) 

and is identical with the two-point Green function of free fields: 

dSp (^1, ^2) = 4l) E n (- |e-^^o+/5..A..y^^(o) (^^) ^(0) (^^)| (5 20) 

LOa is the set of all the field operators: 

c.„ = {A;i,7^'^'^a,,^,, . (5.21) 

X is the set of the coordinates: x = (r, x). 

Using the equations of motion and equal-time commutation relations for uja^ , one can easily verify that 
D^p has the form 

^"pD^, = -Scy ■ (5.22) 
See the definition of A^^ in Appendix C. Thus, the quantity 



M(^-) = ^ (5) (5.23) 
^0 

has a typical structure of a generating functional in quantum field theory; namely, its calculation can be 
carried out by the same diagrammatic technique, and it obeys the same functional equations satisfied by a 
generating functional of quantum field theory. In particular, a literal repetition of the reasonings that are 
presented in quantum field theory shows ^9^ that Af obeys the relation 

= gi5(0)Splna„,^^^.) ^ (5 24) 



where M^''' is computed by the "Wick" rules 



M^^') = exp (Li„t) exp ( ^J^D^f^j" ] exp (^, J,) . (5.25) 
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In (|5.25p . the source J" is introduced only for the fields (that is, J° = Ji = =6 =0), while the Green 
function of the fields satisfies the relation 

K^^pD^p^ = -6^^ (5.26) 

(the definition of A^^ is given in Appendix C). Besides, in the course of taking a variation of gauge 
parameters, M^^*' changes as follows: 

SM^^'^ = -e-^'"' J^^ dr J d^xSHexp Qj^^'f^jJ'^^ e«'^* = (5.27) 
= --^ /^e-'^"°+^^''^'Tr J^^ dr J d^xSHexpiQ)^ , (5.28) 

5M^'^ = e^'°' J^^ dT J d^x 6L exp Q J"!)^^ j/^^ e^'^' (5.29) 

(the sources in (j5.29p are introduced only for the fields) . Relations (|5.28p and (|5.29p are exact analogues of 
the corresponding relations (jC.9P and (jC.5P of quantum field theory. The only difference between M^'*'') and 
a generating functional of quantum field theory arises if one defines the propagators D^^'^^ as the integral 
operators A*^^'^)^^, whose unique definition requires that one should impose certain boundary conditions. 
The propagators for arbitrary are computed in Appendix D. 

As has been mentioned in the beginning of this section, it is necessary (in order to ensure the possibility 
of deducing the Ward identities) that the fictitious fields should have only even-valued frequencies. For the 
remaining fields, we assume the usual relation between statistics and the parity of frequencies. As shown 
in Appendix D, for all the particles, except the fictitious ones, Ai must be chosen as follows: 

Xi = 1 for particles with positive metric , 

(5.30) 

Xi = — 1 for particles with negative metric . 

Notice that such a choice of A^ ensures the coincidence of a matrix element (. . .) with the definition of the 
trace of an operator (in the subspace of the mentioned particles), so that the relation between the parity 
of frequencies and the statistics of a field can be found without a manifest calculation of the propagator. 
For the fictitious particles, the parameters A^ must be chosen as follows: 

Xi = — 1 for particles with positive metric , 

(5.31) 

Ai = 1 for particles with negative metric . 

The series expansion of a fictitious field in the powers of creation and annihilation operators is made 
in Appendix A, where it has been shown that the field C" contains two particles. Correspondingly, Ai are 
equal to 

Notice, once again, that the choice of A.; for the fictitious particles in accordance with (|5.30p implies that 
their Green function contains only odd-valued frequencies. This is a consequence of the fact that the choice 
of A; according to (|5.30p implies that a matrix element (. . .) coincides with the trace of an operator, whereas 
in this case the propagator of any Fermionic field contains only odd- valued frequencies. 

Consequently, given the choice of A^ in accordance with (|5.30p . (|5.32p . a literal repetitior|f| of the 
reasonings of Section 3 makes it possible to conclude that M^^''> obeys the following Ward identity: 

aabt' (x) M(^') = -i^ (^^~PHo+p,.mj.^ J ^yQb^ (^) (y) (x) exp (Q)^ = 

= « / dyQ'ji (y) D'^ (y, x) M^^') , (5.33) 



•^Within the temperature techniques, there remain vaUd the commutation relations and definitions of canonical momenta 
| |2. 171 1- 112. 271 1 under the replacement dt — > idr ■ This replacement also implies the coincidence of the equations of motion 1 12. 291 1 
and 1)5. 9| l. with obvious modifications due to the presence of the operators A'^^. 
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where / dy denotes 

J dy = J^^Ty J d^y . (5.34) 

Besides, one ought to remember that in terms of the variables t the /^-function is defined by the equation 

T'^' (x) D"' {x, y) = i8 (r, - Ty) 5{x~y) (5.35) 

and there holds a relation of the form (I3.17p . 

A variation of the gauge parameters changes M*^^'^ as follows: 

<5M(^') = i [ dxdyQ% (x) D''^ {x, y) li (y) M^^''> , (5.36) 



where A is given by (|4.9p . 

Thus, the statistical average of gauge-invariant operators, and, in particular, the partition function, are 
gauge-independent on condition that 2^q'^''' should be gauge- invariant. 
To prove the gauge invariance of Zq^' ' , let us apply the formula 

5e->^Tr = -e-^TJ f dr e^"SHe'^^ ^ e'^^T, f dr SH (r) , (5.37) 
JO Jo 

H = Ho~ fi.,N, . (5.38) 



Thus, the variation of IuZq'*"-' equals to the variation of M*^'^*' for the free theory (let us denote the 

(5.39) 



latter by Mq^^') with the vanishing sources: 



j°=o 

From (|5.36p it follows that IhZq^^^ is gauge- invariant. The gauge invariance of InZp'*'''' can also be estab- 
lished by a direct calculation using the expression for Hq obtained in Appendices A and E. 

In physical gauges (the Coulomb gauge x'^^ = — 0, k"^^ = 6ik, aab oo; the unitary gauge 
xf oo, that corresponds to the gauge B,n = 0, with B„i — Q^'fi, being Goldstone Bosons; see, e.g., [11]), 
non-physical particles are absent, so that the proposed definition of the partition function and statistical 
average, given by 

Z = Y.Y[xr{n\e-^''^'^'^^^^i..-)\-), (5.40) 

n i 

where \n) and are defined in accordance with (|5.16p . (|5.32p . while (...) stands for the T^-product of 
operators whose statistical average is to be found, coincides with the usual definition in physical gauges. 
Since the partition function and statistical average of a gauge- invariant operator are gauge- independent, 
definition (|5.40p yields a correct result in any gauge. 

Acknowledgement The author is grateful to A. Linde and E.S. Fradkin for useful discussions. 

A Appendix 

We are now going to deduce the expressions for the canonical momenta conjugate to C and C+, as well as 
for their commutation relations with the help of Schwinger's action principle ^6,. Let us present the part 
of the action that corresponds to the fictitious fields: 

t2 



I dx^-d^C+^K^^Vf + igC+-'K'iT\j^jC'') . (A.l) 

We now find a variation 5W: 



ti 



W = r dx Uc+''T''^C^ + C+''T''^T\jifij5C^'\ + G 

J ti 



(A.2) 
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where 



G 



(A.3) 



In accordance with Schwinger's action principle, G is the generator of variations of the field variables at a 
fixed moment of time: 

SC" = i [G, C] , - i [G, C+'^j . (A.4) 

Then, introducing notation (|2.25p . we deduce from (|A.3P and ()A.4[1 that the anticommutation relations 
p.27p hold true, whereas the anticommutators between He and Ilc+ , G and C+ are equal to zero. 

Let us now prove the fact that definition (|2.25p is, in a certain sense, unique. Namely, we define the 
canonical momenta as 

n^. = -a>f"^(9^C+'^ , n^.^ = /3>f°^9^C'' (A.5) 
and suppose that p.27p holds true. We then construct the Hamiltonian 



xoo^f] a a(3 



1)0 



(A.6) 



Let us now demand that the Hamiltonian equations of motion in the form p.29p for the fields and momenta 
should yield equations (|2.33p and definitions IjA.sp . We have 



C = 



1 



(3 a a(3 
By comparison with (jA.Sp . we find 



— H \ ~ ai [gx A„ + X G 



ai = 1 , /3i = - . 

P 



Next, 



whence 



1 



00 ^^c 



00 ' 



dkG+^ , 



/3 = 1 = /3i , (3. 



(A.7) 

(A.8) 

(A.9) 
(A.IO) 



Therefore, all the parameters are defined uniquely, and the expressions for 11 and H coincide with the 
corresponding expressions of Section 2. Let us now verify that the Lagrangian equations are also fulfilled. 
Indeed, 



He 
tie 



1 



1 

^00 



dkG+' + igG+''x\T%^, 



x°''dkU^c - dk ( 



^Ok^O'i. 



1)0 



Vab/^b ■ „ , ^aj^b , „ /^b 



(A.ll) 
(A.12) 



Substituting expressions (|2.25p into (jA.lip . (jA.12p . we find the equations of motion (|2.33p . 

Let us now expand the operator G in the powers of the creation and annihilation operators in the free 
case: 



Lo = -a^C+"x^''9,C+" + G+'^ruabG'' 

mat = -gKit ■ 



(A.13) 
(A.14) 



For the sake of simplicity, we suppose that there exists a matrix S such that transforms the matrix m to a 
diagonal form: 

SabmbcScd = ^J'a5ad ■ (A. 15) 



Introducing the fields 



V+ = G+'-Sba , Va = S-^G"- 



(A.16) 
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we bring Lagrangian (|A.13[) . as well as the equations of motion and commutation relations, to the form 



while the remaining anticommutators are equal to zero. 

From (jA.lSP it follows that the expansion of the field Va has the from 



(A.17) 
(A.18) 
(A.19) 



Vaix) 



dp 



[e-^^'^da ip) 6 {po - <) + 5+ (p) e'P^d {po - n;)] , (A.20) 



1 m^^h' 



00 



,0i , -Ofe 



- 1 p^pk + 



Ma 



1/2 



^00 I ^0 



I PtPk + 



Ma 
^0 



-,1/2 



(A.21) 



Of course, >f'^'^ and /Za must be subject to such equations that and il^ should be real-valued. 
With respect to we assume 

- (^): • (A.22) 
In the general case, C"*" is not Hermitian-conjugate to C. Then, the canonical commutation relations lead 
to the following rules for the operators d and C: 



-{da {p),d+ iq)}^{ba {p),d+ iq)} 
Let us verify, for instance, (|A.19p . namely, 



rSabS [p- q) 



(A.23) 



2/0 2 



001 



x'^'^Pu dk dp 



(2^)^ [^^'^p^^^.p^^ 



1/2 



-ipx-\-iky 



\Zm\ 



5 (po ~ n;) 6 (fco - nt) - 



(po - u;;) S (fco - ujI) 5 {p- fc) = iSabS (f - y) 



(A.24) 



The expansion of the initial fields has the form 



C (x) 



1 



dp 



(27r)' 



—Sac (po - 4 (p) + 



+ e'p-s {po - n;) bt (p)] , 



dp 



{2TTf -^xO-p, 



[e^P^S{po-u;^) dtip) + 



+e-^P^6{po-n;) bc{p)]S-:. 



(A.25) 



(A.26) 



Relation (|A.23[) shows that the field C contains two kinds of Fermions; one of them has a positive norm, 
while the other one has an indefinite norm. Of course, this fact is in agreement with the theorem on the 
relation between spin and statistics. 

Let us, finally, present an expression for the Hamiltonian in terms of the creation and annihilation 
operators: 



Ho -- 
= fd'p 



d^3 



djC + C+'^mabC'' 



.00 



— {p) da ip) + ^^^;b+ ip) ba (p) 

I \^ I 



(A.27) 
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The propagator of the field C is given by 

(0 \TC'^ (x) {y)\ 0) = Sac (0 [TV, (x) V+ (y)] O) 5,^ 



'^P ^-^p{^-v)Qab ^ (A 28) 



(2-)^ 



G'S (P) = ^ad— S-^^ = {x^^y^p.Sab - mabV^ . (A.29) 

B Appendix 

Let us now prove that the second term in the r.h.s. of p.l2p equals to zero. 

We use the antisymmetry property for a product of the fields C° at coincident points, as well as the 
possibility of freely integrating by parts, since the additional terms arising due to the non-commutativity of 
do and T-product are proportional to /'"^■^Sm, which is equal to zero. We then have the following equalities 
(we imply integration over dy yet do not present it explicitly): 

_g2j^dcbj-fbnj^c^^^uQ^(j+f^n(jd ^ 

J ^ jdch jhfn _ jnch jbfd'^ jj^M^^ (J~^^ C^C^ 

= -y C+' d ^>c^''' f'^^Alf'^fC'^Cf . (B.2) 
Here, we have used the Jacobi identity 

jdcb jbfn _j_ jncb jbdf jfcbjbnd 

Next, 

^2 



-,2 



= i^-C+'>^^T%^,f'fC'Cf . (B.4) 

In (|B.4[) . we have also used (12. 4p . Summarizing (|B.1[) . (jB.2[) and (|B.4[) . we find that the second term in the 
r.h.s. of p.l2p has the form 

^gC+^T^^f'^fC'^Cf, (B.5) 
which equals to zero owing to the equations ()2.33p for C+ . 

C Appendix 

Let us now deduce a formula for the variation of Z\\r corresponding to a variation of the parameters in the 
Lagrangian. In the interaction representation, Zw can be written as follows: 

Zw = (0 \Tw exp {iQ + iLmt) \ 0) . (C.l) 

When taking a variation of the parameters in the Lagrangian, we need to examine a variation of the 
vertices and propagators. According to (jC.ip . the variation of vertices is given by an insertion into the 
Green function of a "field" i J dx 5 Lint ■ 
Since the field propagator equals to 

where A-^ is a differential operator in the free Lagrangian, 

Lo = , (C.3) 
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its variation reads 

SDf^ = Ai. i^SAj^i) Dt, ■ (C.4) 

From (jC.4p it follows that the variation of propagators is equivalent to an insertion of the "field" i J dx SLq. 
Thus, the total variation of Green's functions with respect to the parameters of the Lagrangian is given by 
an insertion of the "field" i J dxdL, namely, 



w 



Tw J dxSL^ exp {iQ + iLi, 



If we do not pass to Wick's rules, Z is given by the formula 

SZ^{0\Twexp iiQ~iHint)\0) 



(C.5) 



(C.6) 



where T denotes the symbol of the usual T-ordering. Let us introduce a unifying field, a;^ = (tt^, ipi). Then 
(|C.6p is computed by the usual Feynman rules with vertices defined by —iH-mt and by the propagator 



with being a differential operator in the expression 

1 



(C.7) 



A" 



A(i) 



(C.8) 



where A^^^ is determined by iTiipi and is independent of any parameters; A*-^-* is determined by the free 
Hamiltonian. Next, a literal repetition of the reasonings that have lead us to formula (jC.Sp yields 



5Z= (Q 



T^i j dx SH^ exp {iQ - Hint) 



(C.9) 



We also note that a variation of the Lagrangian and that of the Hamiltonian with respect to the parameters 
are related by 

- SHl,^ = SL\.^^ , (C.IO) 

where the notation 

stands for a variation of A with a fixed u. 

Formulas (jC.SP and (|C.9p have the same appearance. One must, however, bear in mind that the time 
derivatives in (|C.6p commute with the symbol of TvK'Product, whereas in (jC.9p the time derivative that 
arises after the substitution tt, ~ ipi does not commute with the symbol of T-product. Of course, the final 
results, calculated by formulas (|C.5[) and (jC.9|l . are identical, with allowance for relation 



D Appendix 



We are now going to compute the partition functions and Green's functions, as well as to prove Wick's 
theorem for various free fields with a generalized definition of statistical average. 

Let us examine, first of all, the case of one degree of freedom. 

A) Bose system: 



H = Hq — jiN = aLija~^a , [a"*", a] = a , ~ 1 . 
a(r) =e-"^a, a+ (t) = e"^a+ , 



We now define a generating functional. 



e"^^T,expi / dT[j{T)a+{T) + J+ir)a{T)] 



where 



(D.l) 
(D.2) 



(D.3) 



(D.4) 
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Let us transform the r7.-product in (|D.3|) to the normal product: 

exp 1^'^ rfr (J (r) a+ (r) + J+ (r) a (r)) | = cxp dr^ dr^J^ (n) A (n, T2) J (r2)| 

X : exp { Jia+ + Ji^a} : , 



A(Ti,r2) = (ti - T2)e 



-cj(ri-r2) 



7-+ „-i^T 



Ji = / dr J (r) e"^ , Jr = / dr J+ (r) e 
Jo Jo 

Substituting (|D.5P into (|D.3p . we obtairEI 



(D.5) 
(D.6) 

(D.7) 



00 00 -.^ 



k T+k 



,J+AJ 



n=0 k=0 



(aJiJ+) 



n\k\k\ 



A:!fc! 

k—0 n—k 

(aJiJi+)'' d 



E 



X 



k E ■ 



klkl dx^ 

k=0 n=0 

^ kl klkl 

k=0 



c—aXe 



1-X 



= exp I J^^ dTi dT2 J+ (ri )D{ti,T2) J+ (t2) I , 

D (n, = aO (n - .2) e--(--) + - ^ 



Next, the partition function Zq^' reads 



1 



1 - aAe-"^ 



(D.8) 

^e-.e„(n-r.)^^(^^) ^ (D.9) 
n 

(D.IO) 
(D.ll) 



As a result, the generating functional takes the form 



zW(J) = Z«expQj+I?ijj , 



where Di is given by (|D.9p . (|D.10p and equals to 

1 



D {TUT2) = ^ 5] A" (n e-f^'Tra (n) a+ (t^) 



.(1) 

^0 n 



(D.12) 



(D.13) 



Consequently, we can see that the proposed extension of statistical average obeys Wick's theorem (which 
is also true for arbitrary A). Besides, the choice 



X = a 



(D.14) 



implies that the Di-function contains only even frequencies. As has been observed in Section 5, in this case 
there holds the relation 

^a"(n|Qig2|n) = ^ a" (n IQ2Q1I , (D.15) 



being a consequence of the identity 



l = J2\n) a" {n\ 



(D.16) 



*We remind that £n = 7rn//3. 
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which makes it possible to find the spectrum of frequencies of the function Di by using the Bose-properties 
of the operators. Nevertheless, the choice A = —a implies that the function Di (the Green function of Bose 
operators) contains only odd frequencies. It is interesting that the distribution function in this case is also 
of the Fermi character: 



N=^ Y. \e-^"aa+a\ n) = {f^ + l)"' . 



7(1) 



(D.17) 



In a similar way, we deduce the expression for the generating functional of Green's functions for the field 
(r) = {e-'^'^a + e'^'^a+), namely, 



-^"Tr exp (^J^ dr Jip^ = Z^^^ exp Q J£)2 , 



(D.18) 



D2{n,T2) = e-f'''Tr<p{n)<p{T2) u) = - ^ fi-^-^^^-^^'l^a (£„) , (D.19) 



.(1) 

■^0 n 



B) Fermi case: 



H = aujo^a , {a, a"'"} = a, c? = 1, 

Tr exp I J'^ dr (77+ (r) a (t) + a+ (r) rj (r)) | = e''^'^'' : exp (77+a + a+771) 

vt = dre-"V(T) , Vi= dTe'^''r]{T) . 
Jo Jo 



(D.20) 

(D.21) 
(D.22) 

(D.23) 



Calculating the generating functional, one should take into account the anticommutation properties of a 
and rj: 



We have 



n e 



-(3H 



1 + aXe-'^'^ 



(D.24) 
(D.25) 



-^"Tr exp (^j^ dr {r]+a + a+r;)^ 



1 + aXe-'^P 



= e"""^" [1 + aXe-/^"^ - Xrttme~'^^ (l |a+a| l)] = exp {ri+Dsri) 
Da (n, T2) = ae (n - T2) e-'^f^^-^^^ ' 

a 1 — aA (— )" 



£'3 (En) = - 



2 ze„ — w 



(D.26) 

(D.27) 

(D.28) 



We can see, once again, that the choice A = a implies that the Green function, just as it should be in the 
case of a Fermi Green function, contains only odd frequencies. However, in the case A = —a a Fermi Green 
function contains only even frequencies. This definition leads to a Bose distribution function of a Fermion: 



(2) 
^0 n 



(D.29) 



Wick's theorem, once again, is valid for an arbitrary A. A generalization of the above reasoning to the case 
of a system of particles in a space of arbitrary dimension is evident. 
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We finally compute the partition function and Green function of a fictitious particle, whose field operator 
is given by the formulas [see, ((X25)) . (|X26)) and (|X27)) ] 



C"' (x) = SabVb {x) , (x) = V+ (x) Sba , 



Va (X) = 
Va^{x) = 



dP'p 



1 



ik 



->C ] P^Pk + ^^la 



We restrict ourselves to the following values of \d and Afe: 



-''v^-'P^da (p) + e">+'P^b+ (p) 

e'^>+'P^d+ (p) + e-^^>-'P^ba {p) 
1 



Af, = —\d — A- 



A^ = 1 



We find 













e^^^Tr exp 


/ dr{i 


l+^C [x) + C+" (y) 0'') 






Jo 



- 2/) = ^ E ^rf'^r |e-''^r.C'' (x) (y) I n) 



(D.30) 
(D.31) 

(D.32) 
(D.33) 

(D.34) 

(D.35) 
(D.36) 



Zo 
d^p 

!-(-)" A 

2 {}i'^''Pf,p,y - fla) 



(D.37) 
(D.38) 



In (jD.38[) . one needs to make a replacement: po iSn- We can see that the choice A = +1 leads to a Green 
function (|D.38p that contains only odd frequencies. However, the case A = — 1, corresponding to the choice 
of \d and Xb indicated in Section 5 [see, (|5.3ip . (|5.32p ]. leads to a Green function of Fermionic fictitious 
particles that contains only even frequencies. 

Using the reasonings that have been presented in this appendix, one can easily see that the generating 
functional for an arbitrary field linear in the creation and annihilation operators can be computed with the 
help of Wick's theorem: 



n) — Zne'^ 



\JDJ 



where D is the Green function of the field uj: 



(D.39) 



(D.40) 



This makes it possible to prove (|5.17p and (|5.20p . In order to prove (|5.22p . one has to use the equations of 
motion and canonical commutation relations for lo. 



E Appendix 

Let us now perform a canonical quantization of the system of free fields A^, being Goldstone Bosons. In 
view of a tedious character of the resulting formulas, we restrict ourselves to the case >cf £,f , aab = ctSab, 
>c^^ — g^^ . Then the system which consists of a multiplet of vector and scalar fields is equivalent to a set 
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of Abelian vector fields, each interacting only with its scalar field. The Lagrangian of this system of fields 
il/j = {Afj^, cr) has the form 



L 



4 A"' 



Af2 



Let us present the field as follows: 



(E.l) 



A,. 



(E.2) 



Substituting expansion (|E.2p into the equations of motion for and a, that follow from (jE.ip . we find 



□ + m2) v;, - , 

n-—jip+i/3M+ — ja=il + ^]Dip^ 



From (IE. 41) it follows that 



7\f 



(□ + /3M)^ (y5 = (□ + l3Mf CT = . 



(□-a/32)CT = 0. 



(E.3) 
(E.4) 

(E.5) 

Consequently, we can see that is a usual massive vector particle, whereas the fields (p and ct have 
to be decomposed into summands whose Fourier transformations are proportional to S (fc^ — /SM) and 
S' {k"^ — j3M^. Substituting this decomposition into (|E.4[) . we can see that amongst the 8 amplitudes there 
are only 4 arbitrary ones. In addition, the fields and ct can be presented in the form 

^ = / ^^J0^O (fco) {e-^'-" [5 (fc2 - PM) ck + {(3M + ^) S' {P - pM) (cfe - 4)] + h.c.} , (E.6) 

^ = / ^^0^0 (fco) {e-'"^-" [6 ~ PM) dk + [f3M + ^) ,5' - (3M) (c,. - 4)] + h.c.} . (E.7) 

In quantum field theory, Cfc, dk and c^, are operators. Their commutators can be found due to the 
commutation relations between the fields and conjugate momenta. Let us present a qualitative analysis of 
finding these relations. The fact that ct and Her — a must commute with A^ implies that ct and ct commute 
with ip, whence it follows that [c, c?+] — 0. In order that the commutator function should not contain 5", 
it is necessary that the equality [c, c'^] + [d, d'^] = must take place. Finally, the canonical commutation 
relations of ct and Her yield [rf, d'^] = 1. 

Consequently, we find that the field flp is presented in the form (jE.2p . (jE.6p . (|E.7p . where 



d*x\/2ko ^ 
(27r)^/^ 



e (fco) {e~'''^S - M^) 44 + h.c.} , 



(E., 



u\. being three orthonormalized transversal polarization vectors. The creation and annihilation operators 
obey the commutation relations 



7' n+'-' 



Sii'S (k ~pj , [dfejdp] = - [c/c,c+] = S (^k ~p 



(E.9) 



the remaining commutators being equal to zero. By direct calculation, we can prove that all the canonical 
commutators between the fields flfj and the canonical momenta constructed from Lagrangian (jE.ip are 
fulfilled. 

The Hamiltonian of the system equals to 



H = J d^k 
f3M 



2nk 



dt) {ck ~ dk) 



i^k 



fc2 + [3M , 



(E.IO) 
(E.ll) 



The (5'-function can be presented as -^d/dkoS {k^ — /?M) and then integrated by parts. As a result, we 
find that the fields ip and ct can be presented in the form (|E.6p . (jE.7p . with the following replacement: 



.'(.^-,M)^(^ + |-).(.^-,M). 



(E.12) 
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This representation is useful in the case /3 = 0, when the J'-function needs an additional determination. 
In the general case, described by the Lagrangian 



4 



l\f2 

_0_^a,a.,^ 

2 

h 



L 

r = x^'-a^Al^ + w^ (E.13) 
one should proceed in a similar way. Let us now introduce a field V^, 

= - Tl-^M^" . = ^^^^a^K" + {^^^^d^dJab + Ma^ab) a' . (E.14) 

In terms of this field, the equations of motion acquire the form 

Altv''''' = , Ard^V^''' + xbaabct' = , 

A^t - Sab [g,. (□ + m!) " d^d,] + d,d,>cj , d, = . (E.15) 
We now have to find the spectrum of the system. For this purpose, we need to calculate det A, namely, 
In det A = Sp InA = Sp InA"^ (1 + ^A) I 



Sp InA + X! ^'^ 

k 

dud I, 



a 



(□ + Af2) Sab 



^MXb + B^dAIaKb 



Sp InA + -Sp 5^,5ab In [1 + A (□ + Ml] 



1. 



-Sp(5„^ In 



1 + \Mx 



xM 



X — OQ 



= det -^xab det (df.d'' + XabMb) [| Ma det ^ (□ + M^) . 

a 

Thus, the spectrum of the system is determined by solutions of the equations 

□ + Ml = , + ^afcMfc = , 



(E.16) 



(E.17) 



where ^a stand for the eigenvalues of the matrix iriab — MaXab- Taking account of the expression (|E.14p 
for t°, we can see that the field has to be decomposed in 5 (/c^ — M^) and 5 {x^^'^k^ki, — Ha), whereas 
the field a should be decomposed in S {k^ — Af^), 5 {x^'^kf^k^, — fia) and 5' {x^^'^k^k^ — fia). 

Let us also present the expressions for (Wick's) field propagators. The simplest way to find them is to 
calculate the Gaussian functional integral of L. In addition, we can see that if one introduces sources for 
the fields and i° (rather than those for A" and cr") the integrals over and t°- factor out and can be 
easily computed. As a result, we have 



(F;r) = o, 

(i" (x) {v))=ta'^H[x-v) . 
Hence, the propagators of the fields and cr" are found with the help of the relations 
a- = {x^^'^d^^d^Sab + mabV^ Mb {t' - xf^^d^V,') , 



(E.18) 

(E.19) 
(E.20) 



= + TT^t^ (^''"d^dJab + niaby' Mt {t" - x'^'^dX) ■ 



(E.21) 



F Appendix 

In this appendix, we deduce, for the sake of completeness, the Ward identities in case the Green functions 
include the fictitious fields. Deriving these identities is quite similar to the corresponding calculation of 
Section 3, and, therefore, we are not going to present a detailed analysis. 



20 



Consider the generating functionals 



Texp 



Texp 



where 



Q = Q + e+^c 

An analogue of formula p.l2p now takes the form 



C 



-\-ana 



(F.l) 
(F.2) 

(F.3) 



dy tf (y) af,T^''Zcc+ ^ dy { [39+" (y) /"'^"C" (y) - Qj, [y]] {y) (y)) + 



g>c^-\7f (d.C+' (y) f"' ^ {y)) + 9^0+" {y) K^VfT'^^ {y) {y)] C {y) {x) 



dy 

In this appendix, we use the following notation: 

(...) = (r(...)exp (iQ) 



(F.4) 



(F.5) 



As in Section 3, we need to carry out an integration by parts in the second term of (jF.4[) . We, however, 
must take into account that Q depends on C and , and, therefore, does not commute with the symbol 
of T-product. We need to extract from the symbol of T-product, then integrate by parts, using the 
antisymmetry of the operators C (as has been done in Appendix B; one only has to remember that now 
stands aside from the symbol of T-product), then integrate by parts once again, and finally insert da in the 
symbol of T-product. In transposing with the symbol of T-product, we need formula p.3p . while also 
making the replacement Q —i- Q. As a result, the Ward identities have a surprisingly simple form: 



laabt' {x)Z = -(^l dyQ\ {y) {y) (^)) + (/ ^yt'' {y) atdO" (y) (x) 
dye+'{y)f''fC''{y)Cf{y)C+^x)) . 



(F.6) 



This identity has been verified by an explicit calculation in quantum electrodynamics (within a gauge 
in which particles are free) as well as in the first perturbative order of a non-Abelian theory. 

It should be noted that the process of deriving the Ward identities with the help of the usual procedure 
of a non-local change of variables in the functional integral TT] for the generating functional leads to a 
much more involved expression than (jF.6[) . It is interesting, however, that there exists another change of 
variables {supertransformation), that allows one to obtain (jF.6|l . as well as the usual Ward identity. 

Consider the generating functional 

Z = j dAdipd^di^dC dC+e'^+'^ , (F.7) 

where the expression for L is given by formula (|2.ip . Let us subject (|F.7p to the change of variables 
(|2.13p - (|2.14p . and let us choose the gauge parameter as follows: 

A'^ [x) = ^iC [x) , (F.8) 

where /i is an anticommuting object. Notice that the choice (|F.8|1 leads to an exact form of transformations 
(|2.13|) - (|2.14p . since /i^ = 0. Besides, let us make a change of the field C^, namely, 



The corresponding variation of the Lagrangian reads 

5L = 



dxC+''6T''''C'' = I dxC+''T''''f'"^fC'^C^ 



(F.9) 



(F.IO) 



where ST""^ is the result of taking a variation of T'^^, whereas a transition from the first equality to the 
second one in (jF.lOp is compensated by the following transformation of the field C: 



C (x) C {x) - If^dbj^d (jb ^ 



(F.ll) 
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where A is given by formula (jF.Sp . As a result, we find that Lagrangian (|2.1[) . or, more exactly, the action, 
is invariant under (super)transformations (|2.13[) - (|2.14p . (jF.lip with the parameter A defined by formula 
(|F.8p . whereas the Jacobian of this change equals to 1. Therefore, variation affects only the term with the 
sources, and so we obtain the following relation: 



Differentiating (jF.f 2p by 5/69'^ {x), we obtain the Ward identity in the form (|F.6[) . 
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